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In [J. Math. Phys. 53, 042105 (2012)], Baym claims to prove the consistency of the purported 
piece-wise solutions to the fractional Schrodinger equation for an infinite square well. However, his 
calculation uses standard contour integral techniques despite the absence of an analytic integrand. 
The correct calculation is presented and supports our earlier work proving that the purported piece- 
wise solutions do not solve the fractional Schrodinger equation for an infinite square well [M. Jeng, 
S.-L.-Y. Xu, E. Hawkins, and J. M. Schwarz, J. Math. Phys. 51, 062102 (2010)]. 



The one-dimensional fractional Schrodinger equation [TJ [5] is given by 

lj{x 



= Z?„ {-teAy^' ^{x, t) + V{x, t) i,{x, t), (1) 



where Da is a constant, A = d"^ jdx^ is the Laplacian, and (— /i^A)"^^ is the quantum Ricsz fractional derivative, 
which is defined as 

1 /'+00 

(-/i^ A)"/2 ^(^^ t) ^ — y e^-^^l^ \vV t) dp , (2) 
where 4>{p, t) is the Fourier transform of the wavefunction, 

/+00 
i^{x,t)e-'P'''^'dx . (3) 
-oo 

When a = 2, the quantum Riesz fractional derivative becomes equivalent to the ordinary Laplacian, and we recover 
the ordinary Schrodinger equation. 

Unless a is an even natural number, the quantum Riesz fractional derivative is a nonlocal operator, and the usual 
piece- wise approach to solving the ordinary Schrodinger equation is inapplicable. Despite this, in Ref. [Tj, Laskin used 
the piece-wise approach for the fractional Schrodinger equation with an infinite square well potential; in particular, 
he claimed that 

, , , fAcos(f^) for Ixl < a 
I otherwise, 

is the ground state solution. In Ref. [3 , we demonstrated that for — 1 < a < 1 and a ^ Q, this '4!{){x) is not a solution. 
If "tpoix) were a solution, then (— /i^A)" ^ V'o(2^) would vanish at a; = a (one of the ends of the infinite square well) 
but we showed that it does not. More specifically, (— /i^A)"^^ i^aia) is proportional to 

/(«) = £ cos^ Qttp^ dp , (5) 

and we showed that df /da > 0, therefore /(a) only vanishes for a = 0. The integral does not need to be explicitly 
evaluated to reach this conclusion. 

In Ref. 4 , Baym claims to evaluate (— /i^A)"^^ -0o(a^) explicitly to (incorrectly) conclude that /(a) = for all a. 
To clarify this matter further, we not only point out several errors in Ref. [4 , but also provide the correct calculation. 
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I. ERRORS 

Baym begins with the fractional Schrodinger equation ([T]) and claims that with an infinite square well potential, this 
is equivalent, after separating variables, to his Eq. (6) — the free, time-independent, fractional Schrodinger equation 
with the "boundary" conditions ip(—a) = ip{a) = 0. 

These are not equivalent. Indeed, these are not even equivalent for the ordinary Schrodinger equation. The 
solutions of the Schrodinger equation with an infinite square well vanish outside [—a, a], whereas the solutions of the 
free Schrodinger equation with these boundary conditions are sine waves that continue outside of [—a, a]. 

The problem is (at least partly) that Baym has ignored the point of our paper [3] . He claims in his conclusions that 
"the solution inside the well is consistent with the outside." Because the fractional Schrodinger equation is nonlocal, 
it is meaningless to speak of solutions inside and outside. The system must be treated as a unified whole. 

More concretely, Baym considers the integral 

i=r ''^'"T^-'^^^^-^-c^,, (6) 

where q = and claims that this is "a singular integral with poles on the real axis at g = ±1". In fact, the 

integrand is continuous and bounded, because the zeroes of the cosine cancel the zeroes of the denominator. 

The most serious error is his claim that ^ can "be evaluated via analytic continuation as a Cauchy principal value 
integral." After breaking the integral ^ into pieces, he claims to evaluate these pieces by closing the contour in the 
upper or lower half plane. This technique requires an integrand that is an analytic function of q over the relevant half 
plane. 

Most of the factors in the integrand are explicitly analytic functions, but the trouble is the factor of \q\°'. This is 
equal to q" for q > and (—(?)" for q < 0. These functions have obvious analytic continuations, but it is impossible 
to join them together, because they are never equal for q ^ 0, therefore this integrand does not have an analytic 
continuation to the upper or lower half plane. 

In Ref. [S], Baym claims that in Ref. [1] he analytically continued the integrand by the "replacement" of \q\" with 

{iqr + (-»g)° . . 

2cos(a7r/2) ' ^ ' 

This expression is ambiguous. It is only equal to for real q if we take the branch cut for the first term in the 
upper half plane, but the branch cut for the second term in the lower half plane; however, that does not define an 
analytic function over either half plane. There are infinitely many possible interpretations of ([t]) — corresponding to 
the possible choices of i" and (— «)" — but none of them reproduces Baym's evaluation of /. 



II. ANALYTIC COMPUTATION 



It is possible to evaluate the integral ^ in closed form, but this is not as elementary as Baym claims. First, note 
that the integral can be trivially rewritten as 



r+oo 

I{x) — 2 / — cos{TTq/2) cos{TTqx/2a) dq . (8) 

^0 q- — I 

This integrand is easily continued to an analytic function of q by taking the branch cut of along the negative real 
axis. With this, the integral Q can be equivalently taken along a contour, C, that follows the positive real axis, 
except for a small clockwise detour around q — 1; see Figure [l] 

The product of cosines can be written as a sum of 4 exponential functions. It is convenient to break the integral 
([s]) up in this way. Now, define 

J(A)^ / -^e^^-^dq, (9) 

so that 

^(-)-M^(f [f + i]) + ^(f [-f-i]) + ^(f [-! + i]) + ^(f ■ 

In order to evaluate J(A), we will need the related integral 

roo j_a 

K{X) EE I ^^e-^'dt (10) 
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q=l 



FIG. 1: The contour, C. 



for A > 0. To compute this, first note that 



K{0)= dt= tan"6id6l = - ia, i + ia) 

Jo ^ + 1 Jo 

_ r(i-ia)r(i + ia) _ 



2r(i) 



2sin(f + fa) 



2 2 ' 



where we have used the Beta function and Euler's reflection formula. 
Then note that 



cosA + tsinA+ ^e-'^ / e^-'+'^'ds - ^e'^ e'^'^~''^'ds = 



Inserting this identity into the integral (10) and reversing the order of integration gives 



poo poo ^a+i 

-ft'(A)=cosA / ^ (it + sin A / ^ dt+i^e 

Jo ^ + 1 Jo ^ + 1 

= f (cos A sec ^ - sin A esc ^) + |r(a + Vje"'^ [ 

Jo 

By changing variables to z = —is, the first integral becomes 



A /'oo 
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"dt ds — ^e" 
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s-"-ie'"ds- fr(a + l)e 
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'e-'^ds 



~a— 1 As 



e'^ds 



-iX 



{-izy-'e-'Hdz = e'2°" [r(-a) - r(-a, -iA)] , 



in terms of the (upper) incomplete Gamma function. 
By angle addition formulae, 



I (cos A sec ^ — sin A esc ^) — tt 



cos A sin ^ — sin A cos ^ 7rsin(A+fa) 



2 cos ^ sin ^ 



smTra 



By Euler's reflection formula. 



T{a + l)r(-a) (^|e~'^e-*i" - ^e^^e'f" 



TTCSC — 7ra 



) sin (A + fa) 



-7rsin(A + |a) 



2 ^ 2 

These terms cancel, leaving 

K{X) = r(a + 1) [|e'(^+f ")r(-a, lA) - f e-'^^+f ")r(-a, -iX) 



smTra 



(11) 



Next, to evaluate J(A), first suppose that A > 0. The integrand in ^ falls off rapidly enough in the upper half 
plane that J(A) can be equivalently computed by integrating along the positive imaginary axis. Setting q = it, this 
gives 

roo 



J(A) 



FIG. 2: The integral I{x) for a — ^ and a = 1. Note that this does not vanish at a; = 1. 



For A < 0, the calculation is very similar, but J(A) equals the integral along the negative imaginary axis, plus the 
result of integrating clockwise around q = 1, which is — vri e^'^. Together, these give 



J(A) 



-ie'i^KiX) A>0 
+ ie~^t"/f(-A) A<0 



-iX 



Now, we can return to the original integral. For x > a, 

= sin^[X(|[f + + 

= i sin r(a + 1) (eT(^+") [r{-a, f [f - 1]) - r(-a, f 



c.c. 



For X < —a, /(x) = /(— x). For —a < a; < a, 

lix) = ^eT(-f-i) + ^eT(-f-i) + i + ^e-'^") [if(f [1 + f ]) + if (f [1 - f ])] 



2 ^ '2 

= -'^cos If + sin ^ [K{^[1 + f ]) + if (f [1 - f ])] 

= -TTCOsff - isin^r(Q! + 1) [e'5'(f+")r(-a, f [f + 1]) + e'^'-- 



r(-«,f [-! + !]) + c.c. 



(12) 



If, '00 (2;), the naive ground state solution to the fractional Schrodinger equation given in Eq. Q and Ref. [T] were 
correct, then (12 1 would be proportional to cos but it is not; see Figure [2j 
In particular, 

/(a) = ii(a) = i sin ^ r(a + 1) (i e"^ [r(-a) - r(-a, tt^)] + c.c.) 

is the function that we proved to be monotonically increasing, and which Baym claims is identically 0. This function 
is plotted in Figure [3j 

Although this computation demonstrates quite explicitly that ^'0(2;) is not a solution of the fractional Schrodinger 
equation for an infinite square well, it unfortunately does not tell us what the correct solution is. That problem 
remains open. 
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FIG. 3: The function f{a 



:), which Baym claims is 0. 



